Janossy Densities II. Pfaffian Ensembles. * 

Alexander Soshnikov^ 



Abstract 



We extend the main result of the companion paper math-ph/02 12063 to the case of the pfaffian 
ensembles. 

1 Introduction and Formulation of Results 

Let us consider a 2n-particle pfaffian ensemble introduced by Rains in j^j: Let (X, A) be a measure 
space, 4>i,4>2,... 4>2n be complex- valued functions on X, and e(x, y) be an antisymmetric kernel such 
that 

p(xi,...,x 2n ) = (1/Z 2n )det(()) j (xk)) jik=ly ..2n Pf(e(xj,x k )) jik=lr „ t 2n (1) 

defines the density of a 2n-dimensional probability distribution on X 2n = X x • • • x X with respect 
to the product measure A® 2n . Ensembles of this form were introduced in [S] and . We recall (see 
e.g. [3]) that the pfaffian of a 2n x 2n antisymmetric matrix A = (aj k ),j, k = 1, . . . , 2n, a,j k = —a k j, 
is defined as pf(A) = Y2 T (~ l) sl5n ^u'i x • • • x ai n j n , where the summation is over all partitions 
of the set {1, . . . , 2m} into disjoint pairs {h,ji}, ■ ■ ■ , {i n ,jn} such that i k < jk, k = 1, . . .n, and 
sign(r) is the sign of the permutation ■ ■ ■ ,i n ,jn)- The normalization constant in Q (usually 

called the partition function) 



'In 



X 



■>t, 



det((f>j(x k ))j,k=l,...2n Pf(e(Xj,Xk))j,k=l,...,2n (2) 



can be shown to be equal (2n)lpf(M), where the 2n x 2n antisymmetric matrix M = (-Mjjfc)i fc=i 2n 
is defined as 



Mjk= (j)j(x)e(x,y)(j) k (y)\(dx)\(dy). (3) 

■IX' 2 



For the pfaffian ensemble (^Q) one can explicitly calculate k-point correlation functions 

p k (xi, ...,x k ):= ((2n)!/(2ra-fc)!) f x2n _ k p(xi, . . .,x k ,x k+ i, . . . ,x 2n )d\(x k+1 ) . ..dX(x 2n ), k = 1, ... ,2n 

and show that they have the pfaffian form (|5|) 

Pk{x\, ...,x k ) =pf{K(x i: x j )) i j =1: ,„ k , (4) 
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where K(x,y) is the antisymmetric matrix kernel 

(5) 



K(r , = ( ZlK^nM^M^Mv) Ei<^< 2 n^(^)^(^)(y) 



provided the matrix M is invertible (by definition (e0)(x) = f x e(x,y)4>(y)X(dy)). If X C M and 
A is absolutely continuous with respect to the Lebesgue measure, then the probabilistic meaning of 
the k-point correlation functions is that of the density of probability to find an eigenvalue in each 
infinitesimal interval around points xi, X2, ■ ■ ■ X\~. In other words 

p k (xi,X2, ■ ■ ■ x k )\(dxx) ■ ■ ■ \{dx k ) = 

Pr { there is a particle in each infinitesimal interval (xj,Xj + dxi)}. 

On the other hand, if a is supported by a discrete set of points, then 

Pk{x\,X2, ■ ■ ■ Xfc)A(xi) • • • A(xfc) = Pr { there is a particle at each of the points Xi, i = 1, . . . , A;}. 

In general, random point processes with the k-point correlation functions of the pfaffian form 
Q are called pfaffian random point processes. (|H1)- Pfaffian point processes include determinantal 

point processes f|lU|) as a particular case when the matrix kernel has the form ( ^ \ where 

K is a scalar kernel and e is an antisymmetric kernel. 

So-called Janossy densities Jkjixx, ■ ■ ■ > %k), k = 0,1,2,..., describe the distribution of the 
eigenvalues in any given interval I. If X C M and A is absolutely continuous with respect to the 
Lebesgue measure then 

<Jk,i{xi, ■ ■ ■ Xk)X(dxi) ■ ■ ■ X(dxk) = Pr { there are exactly k particles in /, 

one in each of the k distinct infinitesimal intervals (xj, Xj + dx{)} 

If A is discrete then 

<Jk,i(%i, ■ ■ ■ x k) = Pr{ there are exactly k particles in /, one at each of the k points Xj, i = 1, . . . , k}. 

See |1] and [2] for details and additional discussion. For pfaffian point processes the Janossy densities 
also have the pfaffian form (see [HI) with an antisymmetric matrix kernel Lj : 

Jk,i(xi,...,x k ) =const(I)pf(L I (xi,Xj))i ij=lj ...k, (6) 

where 

L I = K I (Id+JK I y\ (7) 

J = ( J and const (I) = pf(J — Kj) is the Fredholm pfaffian of the restriction of the 



4 , 

operator K on the interval /, i.e. const(I) = pf(J — Kj) = (pf(J + Lj))^ 1 = (det(Jd + Jx 
Kj)) 1 / 2 = (det(Id — JL/)) -1 / 2 . (we refer the reader to section 8 for the treatment of Fredholm 
pfaffians). 

Let us define three 2n x 2n matrices G 1 , M 1 , M X \ T : 

G] k = [ hix) [ t(x,y)4> k (y)X(dy)X(dx), (8) 



x 
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Mj k = J a <f>j(x)e{x,y)<f> k (y)X{dx)X(dy), (9) 
Mf,} 1 = f <j )j {x)e{x,y)<j> k {y)X{dx)X{dy) (10) 

J{X\I) 2 

(please compare (9)- (10) with the above formula (j3J) for M). Throughout the paper we will assume 
that the matrices M 1 and M x ^ are invertible. 
The main result of this paper is 

Theorem 1.1 The kernel Lj has a form similar to the formula J2J) for K. Namely, Lj is equal to 
Li(x,y) = 

T.i< , k <2n^)(M x \ I )p k {y) T,i< ] , k <2n^)(M x \%(e x \ l( t )k ){y) 
Ei< i)fe <2n( e xv0i)(^)(M x V)-V fe (y) -e(x,y) + £i<^< 2 „(ex\/^0(*)(M x \%' (e xv <t> k )(y) 

where (e x \i4>)(x) = f x ^ e(x,y)4>(y)X(dy). 

Comparing 1)11(1 with (JSJ) one can see that the kernel Lj is constructed in the following way: 1) first 
it is constructed on X \ I by the same recipe used to construct the kernel K on the whole X, 2) it 
is extended then to I (we recall that Li acts on L 2 (I, dX(x)), not on 1?{X \ I, dX(x))). 

This result contains as a special case Theorem 1.1 from the companion paper [3]. The rest of 
the paper is organized as follows. We discuss some interesting special cases of the theorem, namely 
so-called polynomial ensembles {(3 = 1,2 and 4) in section 2. The proof of the theorem is given in 
section 3. 



2 Random Matrix Ensembles with (3 = 1, 2, 4. 

We follow the discussion in (Hj (see also jllj and |12j^. 
Biorthogonal Ensembles . 

Consider the particle space to be the union of two identical measure spaces (V, /i) and (W, fi) : 
X = V U W, V = W. The configuration of 2n particles in X will consist of n particles v \ , . . . , v n 
in V and n particles w\ , . . . , uu n in W in such a way that the configurations of particles in V and 
W are identical ( i.e. Vj = Wj,j = 1, . . . ,n). Let £j, = 1, . . . ,n be some functions on V. We 

define {4>j} and e in Q so that 4>j{y) = 0, v G V, 4>j( w ) = £j( w )i w £ W, j = 1, ...to, 4>j(v) = 
i>j-n-i{v), v G V, 4>j{w) =0, w G W, j = TO + 1, . . . ,2n, and e(vi,v 2 ) = 0, vi,V2 G V, e(wi,w 2 ) = 
0, to i, U72 € W, e(f , to) = — e(io, u) = <5 W , v G V, to G W. The restriction of the measure A on both V 
and W is defined to be equal to /i. Then (^Q) specializes into (see Corollary 1.5. in 9 ) 

p(vi, ...,v n ) = coTOst n det(^(vi)) i j = i > ... jn det(^('Uj))i i:?= i j .., jn . (12) 

Ensembles of the form (|12|) are known as biorthogonal ensembles (see [7j, Q). The statement of 
the Theorem 1.1 in the case (|12|) has been proven in the companion paper j^j. The special case of 
the biorthogonal ensemble ((T2*]l when V = M, £j(x) = ipj(x) = x 3 ~ x , and V = {C | \z\ = l},£j(z) = 
ipj(z) = z 3 ~ , such ensembles are well known in Random Matrix Theory as unitary ensembles, see 
H3 for details. An ensemble of the form (|12l) which is different from random matrix ensembles was 
studied in [Jj. We specifically want to single out the polynomial ensemble with (3 = 2. 
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Polynomial (ft = 2) Ensembles. 

Let X = R or Z, (j)j(x) = xP , j = 1, . . . , 2n, and A(dx) has a density with respect to the 
reference measure on X (Lebesgue measure in the continuous case, counting measure in the discrete 
case). Then the formula (|12[) specializes into 

p(vi, ...,v n ) = const n ] [ (Vi - vj) 2 ] [ (13) 

l<i<j<n 

Polynomial (/3 = 1) Ensembles. 

Let X = R or Z, 4>j(x) = j = 1, . . . , 2n, e(x,y) = ^sgn(y — x) and A(dx) has a density 

lj(x) with respect to the reference measure on X (Lebesgue measure in the continuous case, counting 
measure in the discrete case). Then the formula (0) specializes into the formula for the density of 
the joint distribution of 2n particles in a so-called (3 = 1 polynomial ensemble (see [5], Remark 1): 

p(x 1 , . . . ,x 2n ) = const n Yl \ X i~ X j\ ] [ ^(Xj). (14) 

l<«<J<2n l<J<2n 

In Random Matrix Theory the ensembles (|14j) in the continuous case are known as orthogonal en- 
sembles , see [H]- 

Polynomial (ft = 4) Ensembles. 

Similar to the biorthogonal case (ft = 2) let us consider the particle space to be the union 
of two identical measure spaces (Y, fi), (Z, fi), X = Y U Z, Y = Z, where Y = R or Y = Z. 
The configuration of 2n particles x\, . . . , X2 n , in X will consist of n particles yi, . . . , y n in Y and n 
particles z\, . . . , z n , in Z in such a way that the configurations of particles in Y and Z are identical. 
We define {<f)j} and e so that 4>j(y) = y- 7 , € Y, (j)j(z) = jz^~ l , z £ Z, e(yi,y 2 ) = 0, e^i,;^) = 
0, e(y, z) = — e(z, y) = 5 yz . As above we assume that the measure \i has a density a; with respect to 
the reference measure on Y. Then the formula Q specializes into the formula for the density of the 
joint distribution of n particles in a (3 = 4 polynomial ensemble (see Corollary 1.3. in [S])) 

P(yi,---,yn) = const n f\ (Vi ~ Vjt Yi ( 15 ) 

l<i<j<n 

In Random Matrix Theory the ensembles (|15[) are known as symplectic ensembles, see [H]. 

3 Proof of the Main Result 

Consider matrix kernels 

Ki = —JKj, Ci = -JLj. (16) 

The the relation (JJJ) simplifies into 

d = Ki (Id- Id)- 1 (17) 

which is the same relation that is satisfied by the correlation and Janossy scalar kernels in the 
determinantal case (0, |2j). The consideration of K-i and Ci is motivated by the fact that the 
pfaffians of the 2k x 2k matrices with the antisymmetric matrix kernels Kj and Lj are equal to the 
quaternion determinants ([0]) of 2k x 2k matrices with the kernels /C/, Ci when the latter matrices 
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are viewed as k x k quaternion matrices (i.e. each quaternion entry corresponds to a 2 x 2 block with 
complex entries). It follows from Q and (|16|) that the kernel ICj is given by the formula 

r _ V" M -t( -{^o)®4>k -(e</>j) (e^fc) \ / e\ 
K i~ fc®(e&) J + V J • (18) 



j,k=l,...2n v 

Let us denote by Cj the following kernel 



l<J,fc<2n 



0j ® 0fe 0j (8) (e x \i<f>k) J V 



As above, e0 stands for J x e(x,y)(j)(y). We use the notation 0j ® ^ is a shorthand for (f>j{x)4>k{y)- 
To prove the main result of the paper we will show that £/ = ICj(Id — /C/) _1 (in other words we 
are going to prove that = where Ci is defined in (|T7|) ). The proof relies on Lemmas 1 and 
2 given below. Let us introduce the notation (ei<fi)(x) = j I e(x,y)<p s (y)dX(y). We will show that 

the finite-dimensional subspace 7i = Span < ( e< ^ s J , ( ^ a J , ( Cl ^ s J 1 is invariant 







J s=l,...2n 

under /C/ and £/. The main part of the proof of the theorem is to show that Ci = K,i(Id — /C/) _1 
holds on TC. 

Lemma 1 The operators /C/, Lj leave TL invariant and Cj = JCj(Id — /C/)" 1 holds on Tt. 

Below we give the proof of the lemma. Using the notations introduced above in (8)-(10) one can 
easily calculate 

1C I (f)= £ -((G^M-^J^l) (20) 



j=\,.„1n 



* ( -1 ) - E I 11 ) - ( ) . (21) 



j=l,...2n 



Defining the 2n x 2n matrix T as 



T sk = I <t>s{x) e(x,y)My)dX{y)dX{x) (22) 
Jx\i 



we compute 



(23) 



^( e f)= £ ((^-t)m-%.(^;), 

^ ' j=X,...2n ^ ^ ' 

where (G 1 — T) s k = M^ k = fp 4> s (x)e(x , y)4>k(y) X(dx) X(dy) . One can rewrite the equations (20)-(21) 
as 

K ' ( * ) = E « G ' - (G')')M-% (*.)-( ) , (24) 

«=' ( :t ) = E «tf + <g')')m-'). 3 ( * ) - ( f ) < 25 > 



i=l,...2n 
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We conclude that that the subspace 7i is indeed invariant under /C/ and the matrix of the restriction 

e<f>s \ ( -e(p s \ ( ti4> s 



of K.i on H has the following block structure in the basis 







s=l,...2n 



(G 7 - (G J )*)M _1 (G 7 + (G J )*)M _1 (G 7 — T)M~ l 

j (26) 

-I'd -Id 

(in particular Ran(K,i\u) = Span j ( e< ^ s ) , ( Cl ^ s )\ ). Let us introduce some addi- 



tional notations: 



s=l,...2n 



,4 = (G / -(G / )')M- 1 , (27) 
5 = (G J + (G J )*)M-\ (28) 
C = (G J -r)Af -1 . (29) 

B C 

When a matrix has a block form = I I (as it is in our case) the matrix M x 




(Id — M) 1 has the block form 

(Id- A + C)- 1 -Id (B -C)(Id- A + C)- 1 C(Id- A + Cy 1 

| (30) 

-(Id- A + C)" 1 -Id- (B -C)(Id- A + C)~ l -C(Id- A + C)' 1 

As one can see from the formulas (31)-(33) the invertibility of Id — Ai follows from the invertibility 
of M X \ T which has been assumed throughout the paper. We have 

(Id- A + C)' 1 = M(M +(G I ) t -T)- 1 = M(M X \ 1 )- 1 (31) 
C(Id- A + C)- 1 = (G 1 -T)(M +(G I ) t -T)' 1 =M I (M x \ I y 1 (32) 
(B -C)(Id- A + C)- 1 = ((G 1 )* + T)(M + (G I ) t -T)- 1 = ((G 1 ) 1 + T)(M x \ I y 1 . (33) 

Let us now compute the matrix of the restriction of £/ on H. We have 

7 - V" (M x\ixt ( -{zx\i4>3)®<l>k -(e x \i<l>j) ® (e X \i<f>k) \ ( e x \i \ 

3 ' v **®** <t>3®{zx\i4> k ) ) + {o o J- (34) 

J,/£ — i,...zn 

Similarly to the computations above one can see that TL is invariant under Lj and 



* ( It ) = Ej^-P')')!"^)-')., ( * ) 

- E « r - (G')')(MX\')-% ( e f 

l<j<2n ^ 

^("_t) = j Ej( r +(G')')( 1 ' AV )- 1 w(lt) 

- E « T + (G')')(MA')-i)„ ( £ * 

l<j<2n 



ei4> s 




(35) 







(36) 
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and 

& ( e t ) = E « G/ - T )( MXV r% ( It ) - E - t)(m^)-% ( e ^ 

^ ' j=l,...2n ^ J ' j=l,...2n ^ 

(37) 

Therefore the restriction of £/ to H in the basis < ( e ^, s ) , ( 6 , s ) , ( e/ / f s ) 1 has 

[ V -0s 7 \ -<f>» J V / J s=1| ...2n 

the following block structure 

(T — (G I ) t )(M x \ 1 )^ 1 (T + (G / )')(M X \ / ) _1 (G 7 - T)(M X \ / ) _1 \ 

( 38 ) 

-Id — {T — (G I ) t )(M x \ I y 1 -Id- (T + (G^XM^ 1 )- 1 -{G 1 - T^M^ 1 )' 1 J 

Comparing (|30|) . (31)-(33) and (|38|) we see that Ci = K\{Id — /C/) _1 on TC. Lemma is proven. 

To show that Ci = K,i{Id — /C/) _1 also holds on the complement of TL it is enough to prove 

(wo- 1 "\ ,^ ( ° 



it on the subspaces ( q J > an d ( ^ ) > where Hi = Span(ej4> s )k=i i ... ) 2n and 7^2 = 
Span(^ s )i = i iii ! 2 Ii - The inveribility of the matrix Mj implies that actually it is enough to prove 
Ci = Ki{Id — /C/) _1 on the subspaces ( ^ > an d f [H ) ± J ' ^ ere we use ^ ne standard nota- 

tion (Hi)~ L for the orthogonal complement in L 2 (I) with the standard scalar product (f,g)i = U X 
f(x)g(x)d\(x). We start with the first subspace. 

Lemma 2 T/ie relation Ci = K,i{Id — /C/) _1 holds on ( , . ■ ) . 

The proof is a straightforward check. The notations are slightly simplified when the functions 
{ei4>k, e<pk, k = l,...,2n} are linearly independent in L 2 (I). The degenerate case is left to the 
reader. Consider f s £ (TCi) ± , s = 1, . . . , In such that 

(e^fc, fs)i = (e^fc, <t> s )i, k = I, . . . , 2n. (39) 
We are going to establish the relation for f ^ , which then immediately extends by linearity to 
the linear combinations of f „ | . We write 



*(_/.) - . E -/.)*( -("/' 



j,fc=l,...2n 



j,/c=l,...2n 







"3 



E ( g ' m »i') ( :i ) - ( n ) < 4 »> 

i=l,...2n 



(we have used (|39j) in the second equality) and 



/c, ( ) = .E « G/ - (G ff ) t )^ 1 )^ ( -t ) " ( 6 o /s ) 



(41) 
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Combining (|iU|) and (|4Tj) we get 



Id 



E ((<?' + (G^M- 1 ) 

j=l,...2n 







(42) 



Similarly to (|23|) we compute 



^ ( e f ) = E « G/ - T ) M ~%- ( It ) ( 43 ) 

^ ' .7=1, ...2n ^ 5 ' 

because f I (e I <p s )(x)<p j (x)d\(x) = - fj f s (x)(e I (p j )(x) x 
cher with ()39|) allows us to conclude 
Ki{Id — ICi)^ 1 f *j ^ is almost identical to the calculation of Ki{Id — /C/) _1 



It should be noted that fCi 







d\(x) = for all j = 1, ... , 2n. This together with (|39|) allows us to conclude that the calculation of 

with the only 

see the 



difference that in the former one we have to replace the term 
last equation of (@U). Namely 



by 



eif, 




K^Id-Ki) 



-i 





~fs 



Ki{Id-K!) 



-i 



1 



-fa 



+ 



-fs 



Sj V ~<t>3 



= e (imdA+B^id-A+cr 1 ) 

j=l,-,2n 

- E {l/2){{A + B){Id-A + C)- 1 ) s]X () 

j=l,...,2n 



j=l,...,2n 

where A, C are defined in (27)-(29). At the same time 





7 ^ " 


H 




- sj 


A / 













(44) 







E (^)/*( £ ^)W./.)i 

j,k=l,...,2n ^ J ' 







E (^r% 

j,fc=l,...,2n 



{t4>k,fs)l 















= /C 7 (Id - /C/)- 1 

^ ^ such that {ei4>k, f)i = J I (ej<p} s .)(x)f(x)d\(x) = 0, k,j = l,... 2n. Lemma 2 is proven. 



(45) 



Therefore Ci 



-fs 



, s = l,... In. By linearity result follows for all 



To check (fT7)l on 



™ ) we note that K^Id-K^ ' ' 



o / =Cl \ 



for (7 such that 



J 7 g{x)4>k(x)d\{x) = 0, fc = 1, . . . , 2n, which together with the invertibility of M finishes the proof. 
Theorem is proven. 
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